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Abstract 

On a series of examples from the field of viscoelasticity we demon¬ 
strate that it is possible to attribute physical meaning to initial con¬ 
ditions expressed in terms of Riemann-Liouville fractional derivatives, 
and that it is possible to obtain initial values for such initial conditions 
by appropriate measurements or observations. 


1 Introduction 


Many physical phenomena lead to their description in terms of non integer 
order differential equations. Formulations of non integer order derivatives, 
generally called fractional derivatives, fall into two main classes: Riemann- 
Liouville derivatives and Griinwald-Letnikov derivatives, on one hand, de- 
hned as (Podlubny 1999, Samko et al. 1993) 
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or the Caputo derivative on the other, dehned as (Caputo and Mainardi 
1971) 
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where n — 1 < a < n. 


In this article we deal only with the Riemann-Liouville fractional deriva¬ 
tives. Fractional differential equations in terms of the Riemann-Liouville 
derivatives require initial conditions expressed in terms of initial values of 
fractional derivatives of the unknown function (Podlubny 1999, Samko et 
ah 1993), like, for example, in the following initial value problem (where 
n — 1 < a < n): 


oA"/W + afit) = (t > 0) (3) 

= {k = l,2,...,n). (4) 

On the contrary, initial conditions for the Caputo derivatives are ex¬ 
pressed in terms of initial values of integer order derivatives. It is known 
that for zero initial conditions the Riemann-Liouville, Griinwald-Letnikov 
and Caputo fractional derivatives coincide (Podlubny 1999). This allows a 
numerical solution of initial value problems for differential equations of non 
integer order independently of the chosen dehnition of the fractional deriva¬ 
tive. For this reason, many authors either resort to Caputo derivatives, or 
use the Riemann-Liouville derivatives but avoid the problem of initial values 
of fractional derivatives by treating only the case of zero initial conditions. 

It is frequently stated that the physical meaning of initial conditions 
expressed in terms of fractional derivatives is unclear or even non existent. 
The old and ubiquitous requirement for physical interpretation of such initial 
conditions was most clearly formulated recently by Diethelm et ah (2005): 

“A typical feature of differential equations (both classical and 
fractional) is the need to specify additional conditions in order 
to produce a unique solution. For the case of Caputo FDEs, 
these additional conditions are just the static initial conditions 
..., which are akin to those of classical ODEs, and are there¬ 
fore familiar to us. In contrast, for Riemann-Liouville FDEs, 
these additional conditions constitute certain fractional deriva¬ 
tives (and/or integrals) of the unknown solution at the initial 
point X = 0 ..., which are functions of x. These initial con¬ 
ditions are not physical; furthermore, it is not clear how such 
quantities are to be measured from experiment, say, so that they 
can be appropriately assigned in an analysis.” 
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(Emphasis is ours). This quotation highlights the utmost importance of 
the interpretation of initial conditions in terms of fractional derivatives for 
farther applications in various fields of science. The physical and geometric 
interpretations of operations of fractional integration and differentiation were 
snggested recently by Podlnbny (2002). However, the problem of interpreta¬ 
tion of initial conditions still remained open. 

In this paper we shall show that initial conditions for fractional differential 
eqnations with Riemann-Lionville derivatives expressed in terms of fractional 
derivatives have physical meaning, and that the corresponding qnantities can 
he obtained from measurements. We shall also demonstrate that in many 
instances of practical signihcance zero initial conditions, which are used so 
frequently in practice, appear in a natnral way. 


2 Number of initial conditions, past history 
and memory 

When a physical process can be described in terms of a differential eqnation 
of integer order n, it is well known that n conditions are reqnired to solve the 
system. In this paper only initial conditions are considered. It is also known 
(Podlnbny 1999, Samko et al. 1993) that fractional differential eqnations of 
order a reqnire a* initial conditions, where a* is the lowest integer greater 
than a. This means that if a < 1 as is the case in viscoelasticity when inertial 
effects are negligible, a single initial condition is snfficient. However, one of 
the reasons for the snccess enconntered in describing viscoelasticity by means 
of differential eqnations of non integer order is their ability to describe real 
behavionr, inclnding memory effects snch as are observed in polymers, using 
only a restricted nnmber of material parameters. Such memory effects may 
continue to affect the material response long after the cause has disappeared, 
as observed in stress relaxation after a non monotonous loading programme 
(Heymans and Kitagawa 2004). In such a case a single initial condition would 
appear insufficient to predict material response. 

Here we shall consider only the response of a system starting at t = 0 
from a state of absolnte rest. As a farther simplihcation, we shall consider 
only response to ideal loading programs, snch as step or impulse response. 
The effects of a hnite loading time, of the details of the loading program, and 
of past history will be acconnted for separately in a seqnel. 

It has been shown (Beris and Edwards 1993) that thermodynamically 
valid constitntive eqnations for viscoelasticity are completely eqnivalent to 
analog models containing only elements (springs and dashpots) with positive 
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coefficients. A suitable hierarchical arrangement of springs and dashpots 
gives rise to spring-pot behaviour (described below), either exactly at all 
timescales for an inhnite tree (Heymans and Bauwens 1994), or in the long¬ 
term (or low-frequency) limit for an infinite ladder or inhnite Sierpinski gasket 
(Heymans and Bauwens 1994, Schiessel and Blumen 1993, 1995). In the latter 
case short-term behaviour is similar to a Maxwell model with one element 
replaced by a spring-pot. Therefore the equivalence demonstrated by Beris 
and Edwards can be generalized to models including spring-pots (Heymans 
1996), hence discussion here will be limited to such models. 

3 Spring-pot model 

We shall start with a spring-pot alone, which is a linear viscoelastic element 
whose behaviour is intermediate between that of an elastic element (spring) 
and a viscous element (dashpot). The term “spring-pot” was introduced 
by Koeller (1984), although the concept of an element with intermediate 
properties had been introduced some time earlier. The constitutive equation 
of a spring-pot is: 

a{t) = K oA“e(t) or e{t) = ^ oD^^a{t) (5) 

where a is stress, e is strain and K is the model constant. The spring-pot is 
the viscoelastic version of Westerlund’s “simplest model” (Westerlund 2002). 
If a = 0 the element is linear elastic (Hookean spring) whereas if a = 1 it is 
purely viscous (Newtonian dashpot). Insight can be gained from the response 
of a spring-pot in a few simple cases, using the general relationship (Podlubny 
1999, Samko et ah 1993) 

3.1 Creep or general finite load 

In the case of creep, a stress step ag is applied at initial time t = 0. The 
strain response is hence e{t) = {ag/KT{1 -|- a))t'^. The initial value of the 
strain vanishes, i.e. there is no instantaneous (elastic) strain, only an anelas- 
tic (retarded) response. However, the hrst ordinary derivative of strain is 
unbounded, so that a finite though undehned strain can be reached in an 
arbitrarily small time interval. 

The change of e{t) is described by the fractional differential equation 

„Dfe(t) = ^ (7) 
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In accordance with the theory of fractional differential equations in terms 
of Riemann-Liouville derivatives, an initial condition involving QDf~^e{t) is 
required. This condition can be found by taking the first-order integral of 
the constitutive equation as 


oDr^eit) 


t^o 


oD;'(,a„/K) 


t^O 


In the case under consideration stress is hnite at all times, hence 
oD^^ao ^ Q ~ which leads to zero initial condition for namely 




t^O 


= 0 . 


( 8 ) 


The same considerations apply to a general hnite load o'(f). In the latter 
case the equation to be solved is 


oA“e(«) = (9) 

and the initial condition to be attached to this equation is the zero initial 
condition (jH)). 


3.2 Stress relaxation or general deformation 

The stress response to a strain step Cq is a{t) = (eoi^/r(l — a))t~°‘. The 
initial stress is unbounded rehecting the fact that a spring-pot (just like a 
dashpot) cannot respond immediately to a bounded stress: it has an inhnite 
initial modulus or a vanishing initial compliance. However, relaxation to a 
hnite though undehned stress occurs in an arbitrarily small time interval. 

The change of a(t) is described by the fractional diherential equation 

oDrMt) = Keo. ( 10 ) 

From the known value of Cq we can obtain the initial value (as t approaches 
zero) of Clearly, if QD^°‘a(t) is to be hnite although it is dehned 

over a vanishingly small time interval, a{t) must be unbounded. On the 
contrary, the initial value of oD^°‘a{t) is well dehned and hnite, and that of 
is zero. Thus, contrary to the idea expressed by some authors 
(e.g., Glockle and Nonnenmacher 1991), initial value problems expressed in 
terms of fractional integrals are not better posed than those expressed in 
terms of fractional derivatives. 

If strain increases linearly with time, stress increases as Stress is 

bounded, but the initial values of its integer order derivatives are unbounded. 
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The known strain rate allows us to define the initial value of oDj “cr(t). In 
fact, in this case, zero initial conditions are found both for QD^^a{t) and 

For any general finite strain e{t), following the same reasoning, again zero 
initial conditions are found. 

In all three examples given here, initial conditions expressed in terms 
of fractional derivatives or integrals arise naturally when taking measurable 
quantities into account. 


3.3 Impulse response 


The impulse response is seldom used in viscoelasticity except as a mathemat¬ 
ical convenience, because it is even more problematic to apply a homogeneous 
impulse of stress or strain on a sample than it is to apply a step. However, 
we shall investigate the impulse response following the same reasoning as for 
the step response above. 

Consider an impulse of stress defined as BS(t) applied to the spring-pot 
at time t = 0. After that, the stress remains zero. The strain response is 
e{t) = {B/KT. The initial stress singularity gives rise to a lower- 
order strain singularity, since a spring-pot cannot deform immediately. 

The strain e{t) for f > 0 is the solution to the fractional differential 
equation 

oA“e(«)=0. (11) 


In accordance with the theory of fractional differential equations with 
Riemann-Liouville derivatives, an initial condition involving ^ ^ 

is required. This can be found through integration of the constitutive equa¬ 
tion, as 


Dr^eit) 


t^o 


D;^a{t)/K ^ =B/K, 




which gives the following initial condition to equation (HH): 


Dr^eit) ^ ^ = B/K. 


t^o 


( 12 ) 


In this problem in terms of Riemann-Liouville derivatives B is the initial 
impulse of stress <j(t), e(t) is the strain after application of this impulse, and 
the known impulse of stress yields a non-zero initial condition m involving 
a fractional derivative of strain. This fractional derivative is non zero, well 
defined, and bounded. Note that both strain and its integer-order derivatives 
are unbounded, and its first order integral is zero, so that a meaningful 
initial condition expressing the loading conditions cannot be obtained using 
integral-order derivatives. 
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The physically unrealistic stress response to a prescribed strain impulse 
will not be considered here. In fact, the analytical solution has a strong 
^-(i+“) divergence, reflecting the fact that a strain impulse cannot be applied 
to a spring-pot. 

4 The key: look for inseparable twins 

Now, after introducing the above simple example, let us formulate our gen¬ 
eral approach to interpretation of initial conditions involving the Riemann- 
Liouville fractional derivatives. 

In a general case, when we consider some fractional differential equation 
for, say, U{t), we have to consider also some function V{t), for which some 
dual relation exists between U(t) and V(t). For example, in viscoelasticity 
we have to consider the pair of stress a(t) and strain e(t); in electrical circuits 
the pair of current i(t) and voltage v(t); in heat conduction the pair of the 
temperature difference T(t) and the heat flux q(t)] etc. Functions U(t) and 
V {t) are normally related by some basic physical law for the particular held 
of science. In each scientihc held there are such pairs of functions like the 
aforementioned, which are as inseparable as Siamese twins: the left-hand 
side of the initial condition involves one of them, whereas the evaluation of 
the right-hand side is related to the other. 

This concept is not restricted to the spring-pot treated above, but is 
further applied in the subsequent sections to more elaborate models of vis¬ 
coelastic behaviour. Indeed, a spring-pot is a particularly crude model, which 
has several unrealistic and unphysical characteristics. As pointed out above, 
it has a vanishing initial compliance or an inhnite initial modulus. Viscoelas¬ 
tic solids, on the contrary, have a well-dehned instantaneous modulus. (Note 
that an unbounded initial modulus is no more of a problem when describing 
a viscoelastic fluid than it is when describing Newtonian viscosity: if a step 
strain is applied to a dashpot, the initial stress is also unbounded). At long 
times there is no limit to anelastic strain of a spring-pot: creep continues 
indehnitely. Also, stress relaxes to vanishingly small values. The increase of 
stress in constant strain-rate conditions means that if attempting to describe 
a viscoelastic fluid, steady state flow is never attained. When describing a 
viscoelastic solid, again we hnd an unbounded modulus at t = 0. In spite 
of these limitations, the spring-pot can give an approximate description of 
polymer viscoelasticity in the intermediate time range. Several slightly more 
elaborate models, which alleviate some oversimplihcations of a single spring- 
pot, will be investigated below. 
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5 The fractional order Voigt model 


The fractional Voigt model (a spring and a spring-pot in parallel) is nowa¬ 
days generally nnderstood as a long-term approximation to the fractional 
Zener model and it might seem irrelevant to express concern over initial 
value problems for the Voigt model. However, as the purpose of this note 
is mainly to examine how initial conditions endowed with physical meaning 
may be expressed in systems whose constitutive equations contain fractional 
derivatives, we shall continue to examine the fractional Voigt model. 

The constitutive equation of this model is 

a{t) = Ee{t) + KoD^e{t). (13) 


The Voigt element or associations thereof are considered in viscoelasticity 
modelling to be appropriate to obtain the strain response to a prescribed 
stress program, so we shall investigate only such cases here. 

Assume a stress impulse B6(t) is applied to a Voigt element at time t = 0. 
Then the fractional order equation we need to solve for e{t) {t > 0) is 

Ee{t)+KoD^e{t)=0. (14) 


In agreement with the theory of fractional differential equations in terms 
of Riemann-Liouville derivatives, we need an initial condition, which will 
involve the value of oD^~^e{t) for f —> 0. This condition can be obtained by 
integration of the constitutive equation as 


+ K„Dr'c{t) = oD;'a{t) 


(15) 


The limit of the right hand side is the magnitude B of the stress impulse. 
On physical grounds, the spring-pot cannot deform instantaneously under a 
hnite stress, and, as is the case for a spring-pot alone, any singularity of e{t) 
must be weaker than that of the stress impulse, thus 


oD;h{t) 



0 . 


This can also be found from examination of the behaviour of the left hand 


side of the relationship (HHI): if oA ^e(^) 
is unbounded and equation m cannot 
tion dually takes on the form of 


is non zero, then 
De fulhlled 


oDr"e{t) 
Hence the initial condi- 


A'„Dr'e(i) 


S. 


(16) 


This condition expresses the initial value of the fractional derivative of 
strain, in terms of the stress impulse. We see that the initial 












condition is obtained expressing a fractional derivative of the unknown strain 
in terms of a measurable, physically meaningful value of its “inseparable 
Siamese twin” - the stress. The obtained initial condition is, in fact, identical 
to the initial condition of the spring-pot alone. This reflects the known fact 
that the spring in the Voigt model only affects long-term behaviour. 

Now let us consider creep, i. e. the response to a stress step (Tq applied 
at f = 0. The equation to be solved for the strain e{t) is 


Ee{t) + K — Uo, 


(17) 


and the initial condition for (ED can be found from 


EoD;^e{t) + KoDrHt) = 


where the limit of the right hand side is zero. A bounded stress can produce 
only a bounded strain, so the limit of the hrst-order ordinary integral of 
strain in the left hand side is also zero. Thus the initial condition has the 
form: 


Dr^eit) 


t^o 


= 0 . 


(18) 


Once more, knowledge of a measurable quantity (do) leads to an initial 
condition expressed in terms of a fractional order derivative of the unknown 
(e(f)), its inseparable Siamese twin. 

The case of a general hnite stress program is similar to that of creep. The 
equation to be solved is now 


Ee{t) + K oD^e{t) = a{t), 


(19) 


and the initial condition is identical to the condition m obtained in creep. 


6 The fractional order Maxwell model 

To keep to a simple model while describing realistic behaviour for a viscoelas¬ 
tic solid, a spring expressing instantaneous elasticity must be associated in 
series with the spring-pot. This eliminates the unbounded initial stress in 
description of relaxation. The Maxwell element or associations thereof are 
considered in viscoelasticity modelling to be appropriate to obtain the stress 
response to a prescribed strain program, so we shall investigate only such 
cases here. 

The constitutive equation of the Maxwell model is 

eW = ^ (y{t) + ^ 
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or 


^oO“cr(t) + ^cr(t) = oC’“e(f). 


( 20 ) 


In stress relaxation, a step strain eo is applied at t = 0. Then the equation 
to be solved is 

+ h - (21) 


K 


r(l - a)' 


An initial condition is required, involving the value of 
Integrating the constitutive equation (EOl) and considering t 
proaches zero, we have 


re limit as t ap- 


^ 0 ^ oA = oA 


OL — l 


e{t) 


J 


( 22 ) 


Since strain remains bounded during loading, and a < 1, the right hand 
side inside brackets is bounded and vanishes when t —>■ 0. Since the left hand 
side is a linear combination of positive functions with positive coefficients, 
it can only vanish if each term vanishes. This means that stress remains 
bounded during loading, and hence that we obtain the following initial con¬ 
dition: 




t^o 


= 0 . 


(23) 


Here again we observe that the initial condition on the unknown stress 
arises naturally from its Siamese twin, the known strain. 

Now we shall consider the strain impulse response. A strain impulse of 
magnitude A6{t) is applied at time t=0. Thereafter, the equation to be 
solved is 

^ ^ = 0- (24) 

The required initial condition is obtained as above by integrating the 
constitutive equation I® and considering the limit as t approaches zero: 


1 

E 


oA + 


1 

K 




e{t) 


t^o 


(25) 


The limit of the right hand side of (i2nD is A. Hence the limit of the left 
hand side must also be bounded. This means that the limit of the hrst-order 
integral of stress must be bounded, and the a -f-1 integral must vanish, and 
the initial condition is hnally 


1 

E 


oA (^{t) 


= A. 


t^o 


(26) 
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The singularity in the stress response to a strain impulse is now of the 
same order as that of the strain impulse itself: adding a spring in series with 
the spring-pot has weakened the singularity. 

The strain response to a stress impulse is identical to that of a spring-pot 
alone since the impulse response of the spring vanishes. 


7 The fractional order Zener model 


Among the fractional order models of viscoelasticity considered in this article, 
the most general is the Zener model. Its constitutive equation is 

cr(t) + V oD^ait) = A e(t) + fi oD^e(t), (27) 

where A = Boo is the long-term modulus, fi = K(Eq — Eoo)/Eq, v = [ijE^ 
and Eq is the instantaneous modulus. 

Let us first investigate the response to a stress impulse B5(t) applied to 
the Zener element at time f = 0. Then the fractional differential equation we 
need to solve for e(t) (t > 0) is: 

Ae(f) -l- fioE'^6(t) = 0. (28) 


In accordance with the theory of fractional differential equations, we need 
an initial condition involving the initial value of QD(^~^e(t). Integration of 
the constitutive equation gives: 

oD7^a(t) + uoD°-^a(t) = XoD7^e(t) + fi oD^-^e(t). (29) 

The initial condition can be found by considering equation (Ei as f ^ 0: 


oA Cr(^) + '^oA“ Cr(f) = -^oA ^(f) /i e(^) 


Jt^O 


(30) 


Using considerations similar to those in case of the Voigt model under 
stress impulse, we obtain the initial condition in the form: 


fi oDf e(t) 




= B 


(31) 


As in case of the Voigt model, this condition gives the initial value of the 
fractional derivative of unknown strain, oD(^~^e(t), in terms of its “insepara¬ 
ble twin” - the stress. 

The right and left hand sides of equations EH|l and (EDD are formally 
identical, hence following the same reasoning as above the response to a 
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strain impulse A6{t) applied to the Zener element at time t = 0 will be the 
solution to the equation 


a(t) + z/ oD^a(t) = 0 


(32) 


with the initial condition 



A 


(33) 


This formal equivalence between response to a stress or strain impulse 
reflects the well known fact that the Zener model is the simplest model 
capable of describing response to a stress or strain program equally well. 

In case of creep, i.e. a step-stress a(t) = cxo for cr > 0, we have the 
equation for e(f): 



Ae(f) -|- fioD^e{t) 


(34) 


The initial condition can also be found by considering equation (Ei as 


t 0. 


Following a similar reasoning to that given above for the Maxwell model 
in stress relaxation, we find a zero initial condition to accompany equation 

(P|) : 


= 0 . 


(35) 


This initial condition in terms of fractional derivative of e{t) appeared 


from consideration of its “inseparable twin” 

Similarly, for stress relaxation, e{t) = cq, we obtain the equation for a{t)\ 


a{t) + uoD^a{t) = Acq + fieo — - 

i (^i — a 


(36) 


The initial condition for the unknown stress cr(f). 


„A“-V(()],^„ = o, 


(37) 


appears naturally from consideration of the initial value of strain. 

Let us now consider the case of general load a(t) = cT*(f). The equation 
to be solved for e(t) is 


\e(t) + fioD^eit) = a^(t) + u oD°a^(t) 


(38) 


The corresponding initial condition can be obtained using the following 
procedure. Consider some small t = a. Starting at f = 0, stress a{t) must be 
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recorded until t = a, and based on the recorded values the left hand side of 
the relationship (Ei must be evaluated. The obtained quantity provides an 
approximation of the initial value for the expression in the right hand side 
of (j2ni)- 

In some cases it is possible to find the limit of such approximation as 
a ^ 0. For example, for a physically realisable continuous load we 

obtain a zero initial condition in the form: 


nr^t) 




= 0 . 


(39) 


It is worth mentioning that this procedure amounts, in fact, to the same 
as measuring the initial value of, for example, the first derivative in the case 
of classical differential equations of integer order. From the examples given 
above, it can be seen that for any physically realistic model, zero initial 
conditions will be found for a continuous loading program or even in the case 
of a step discontinuity. Non-zero conditions will only be found in the case of 
an impulse. 


8 Conclusion 

In this note we demonstrated on a series of examples that it is possible to at¬ 
tribute physical meaning to initial conditions expressed in terms of Riemann- 
Liouville fractional derivatives. 

To summarize, expressing initial conditions in terms of fractional deriva¬ 
tives of a function U(t) is not a problem, because it does not require a direct 
experimental evaluation of these fractional derivatives. Instead, one should 
consider its “inseparable twin” V{t) related to U{t) via a basic physical law, 
and measure (or consider) its initial values. 

It is worth noting that the only case where non zero initial conditions 
appeared in our considerations, is the case of impulse response. In other 
cases (including the Zener model under physically realisable load program), 
the initial conditions are zero, and in such cases the use of the Riemann- 
Liouville derivatives, the Griinwald-Letnikov derivatives, and the Caputo 
derivatives is equivalent. 
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